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Summary; 


A  major  goal  of  this  project  is  to  identify  suitable  biomarkers  of  selenium 
chemoprevention  in  human  prostate  cell  models.  During  this  first  year  of  funding  period, 
we  investigated  the  cellular  and  molecular  changes  mediated  by  selenium  in  the  PC-3 
human  prostate  cancer  cell  line.  Please  see  the  attached  article  (Dong  et  ah,  Cancer  Res., 
63,  52-59)  in  Appendix  for  detailed  description  of  the  specific  aspects  of  the  research.  In 
the  original  Statement  of  Work,  we  proposed  to  look  at  the  normal  prostate  cells  and  the 
LNCaP  prostate  cancer  cell  line  in  the  first  year  and  the  PC-3  cell  line  in  the  second  year. 
The  reason  that  we  switched  the  sequence  of  the  work  is  that  PC-3  cells  are  relatively 
easy  to  handle.  We  would  like  to  work  out  the  experimental  conditions  first  with  this  cell 
line  before  we  proceed  to  the  other  two  cell  models. 


Key  Accomplishments: 

>  Selenium-induced  growth  inhibition  in  PC-3  cells  is  achieved  mainly  by  cell  cycle 
blockade  coupled  to  an  induction  of  apoptosis. 

>  A  large  variety  of  potential  selenium-responsive  genes  were  identified  by 
oligonucleotide  array  analysis. 

>  Clustering  analysis  grouped  these  genes  into  early,  intermediate,  and  late 
responsive  clusters. 

>  Since  early  responsive  genes  may  be  important  in  initiating  the  cascade  of  events 
leading  to  the  action  of  selenium,  we  classified  these  genes  according  to  their 
known  functions  in  cell  growth  and  tumorigenesis.  These  include  cell  cycle 
regulators,  apoptosis  controllers,  cell  adhesion  and  invasion  proteins,  signaling 
molecules,  transcription  factors,  and  oncogenes/tumor  suppressors. 

>  The  expression  changes  of  10  genes  critically  involved  in  cell  cycle  regulation 
were  selected  for  verification  by  Western  analysis.  An  agreement  rate  of  70% 
was  obtained  from  these  confirmation  experiments. 

>  From  the  array  data,  we  were  able  to  formulate  a  schematic  diorama  of  signaling 
pathways  that  provide  the  supportive  framework  for  understanding  selenium- 
mediated  cell  cycle  blockade. 

>  The  data  also  provided  valuable  insights  into  novel  biological  effects  of  selenium, 
such  as  inhibition  of  cell  invasion,  initiation  of  DNA  repair,  and  induction  of 
TGF-P  signaling. 
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Reportable  Outcomes: 

>  Publication; 

Dong,  Y.,  Zhang,  H.,  Hawthorn,  L.,  Ganther,  H.E.,  and  Ip,  C.  (2003)  Delineation 
of  the  molecular  basis  for  selenium-induced  growth  arrest  in  human  prostate 
cancer  cells  by  oligonucleotide  array.  Cancer  Res.,  63,  52-59. 

>  Abstract : 

Dong,  Y„  Zhang,  H.,  Hawthorn,  L.,  and  Ip,  C.  (2002)  Delineation  of  the 
Molecular  Basis  for  Selenium-induced  Growth  Arrest  in  Human  Prostate  Cancer 
Cells  by  Oligonucleotide  Array.  Proceedings  of  the  American  Association  for 
Cancer  Research,  43: 167 

>  Presentation: 

93rd  Annual  Meeting  of  American  Association  for  Cancer  Research,  April  2002, 
San  Francisco,  California.  Poster  Presentation,  "Delineation  of  the  Molecular 
Basis  for  Selenium-induced  Growth  Arrest  in  Human  Prostate  Cancer  Cells  by 
Oligonucleotide  Array",  Dong  Y„  Zhang,  H.,  Hawthorn,  L.,  and  Ip  C. 

>  Employment  received: 

Affiliate  Member,  Dept,  of  Cancer  Chemoprevention,  Roswell  Park  Cancer 
Institute,  Buffalo,  NY  14263 
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Delineation  of  the  Molecular  Basis  for  Selenium-induced  Growth  Arrest  in  Human 
Prostate  Cancer  Cells  by  Oligonucleotide  Array1 

Yan  Dong,  Haltao  Zhang,  Lesleyann  Hawthorn,  Howard  E.  Ganther,  and  Clement  Ip2 

Departments  of  Cancer  Prevention  [Y.  D.,  C.  I.],  and  Cancer  Genetics  [H.  Z,  L  H.j,  Roswell  Park  Cancer  Institute,  Buffalo,  New  York  14263,  and  Department  of  Nutritional 
Sciences,  University  of  Wisconsin.  Madison,  Wisconsin  52706  fti  E.  G.j 


ABSTRACT 

Despite  the  growing  interest  In  selenium  intervention  of  prostate  cancer 
in  humans,  scanty  information  is  currently  available  on  the  molecular 
mechanism  of  selenium  action.  Our  past  research  indicated  that  methyl- 
seleninic  acid  (MSA)  is  an  excellent  reagent  for  investigating  the  antican¬ 
cer  effect  of  selenium  in  vitro.  The  present  study  was  designed  to  examine 
the  cellular  and  molecular  effects  of  MSA  in  PC-3  human  prostate  cancer 
cells.  After  exposure  to  physiological  concentrations  of  MSA,  these  cells 
exhibited  a  dose-  and  time-dependent  inhibition  of  growth.  MSA  retarded 
cell  cycle  progression  at  multiple  transition  points  without  changing  the 
proportion  of  cells  in  different  phases  of  the  cell  cycle.  Flow  cytometric 
analysis  of  annexin  V-  and  propidium  iodide-labeled  cells  showed  a 
marked  induction  of  apoptosis  by  MSA.  Array  analysis  with  the  Af- 
fymetrix  human  genome  U95A  chip  was  then  applied  to  profile  the  gene 
expression  changes  that  might  mediate  the  effects  of  selenium.  Gene 
profiling  was  done  in  a  time  course  experiment  (at  12,  24,  36,  and  48  h) 
using  synchronized  cells.  A  large  number  of  potential  selenium-responsive 
genes  with  diverse  biological  functions  were  identified.  These  genes  fell 
into  12  clusters  of  distinct  kinetics  pattern  of  modulation  by  MSA.  The 
expression  changes  of  10  genes  known  to  be  critically  involved  in  cell  cycle 
regulation  were  selected  for  verification  by  Western  analysis  to  determine 
the  reliability  of  the  array  data.  An  agreement  rate  of  70%  was  obtained 
based  on  these  confirmation  experiments.  The  array  data  enabled  us  to 
focus  on  the  role  of  potential  key  genes  ( e.g GADDI 53,  CHK2,  p21IIAFI, 
cyclin  A,C  DK1,  and  DHFR)  that  might  be  targets  of  MSA  in  impeding  cell 
cycle  progression.  The  data  also  provide  valuable  insights  into  novel 
biological  effects  of  selenium,  such  as  inhibition  of  cell  invasion,  DNA 
repair,  and  stimulation  of  transforming  growth  factor  0  signaling.  The 
present  study  demonstrates  the  utility  of  a  genome-wide  analysis  to  elu¬ 
cidate  the  mechanism  of  selenium  chemoprevention. 


INTRODUCTION 

Recruitment  to  the  National  Cancer  Institute-sponsored  SELECT3  trial 
began  in  2001.  This  is  a  Phase  III,  double-blind,  placebo-controlled, 
12-year  trial  designed  to  assess  the  effect  of  selenium  and  vitamin  E, 
either  individually  or  in  combination,  on  the  incidence  of  prostate  cancer. 
The  launching  of  this  trial  is  largely  driven  by  the  milestone  finding  of 
Clark  et  al,  (1)  that  selenized  yeast  supplementation  was  capable  of 
significantly  reducing  the  incidence  of  prostate  (RR  =  0.37),  lung 
(RR  =  0,54),  and  colon  cancers  (RR  =  0.42).  The  SELECT  protocol  also 
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provides  for  the  establishment  of  a  repository  for  prostate  biopsy  tissue, 
blood  cells,  and  plasma.  These  materials  will  be  put  aside  for  research 
discoveries  in  the  future.  One  of  the  secondary  objectives  of  the  trial  is  to 
study  cellular  and  molecular  biomarkers  using  the  banked  samples  and  to 
delineate  their  relevance  with  respect  to  prostate  carcinogenesis  and  drug 
effects  (2).  Despite  the  considerable  public  interest  in  the  potential  benefit 
of  selenium  chemoprevention  of  prostate  cancer,  scanty  information  is 
currently  available  on  the  molecular  targets  or  the  signaling  mechanism 
underlying  the  anticancer  action  of  selenium.  Our  present  study  was 
aimed  at  addressing  this  gap  of  knowledge  with  the  use  of  a  human 
prostate  cancer  cell  line. 

Se-Met  is  the  selenium  compound  used  in  the  SELECT.  It  is,  however, 
not  particularly  suitable  for  mechanism  studies  in  cell  culture.  The  reason 
is  that  Se-Met  needs  to  be  metabolized  primarily  in  the  liver  to  a 
monomethylated  intermediate  for  the  expression  of  its  anticancer  activity 
(3-6),  and  epithelial  tissues  generally  have  a  low  capacity  to  generate  a 
monomethylated  selenium  metabolite  from  Se-Met.  Consequently,  con¬ 
centrations  of  Se-Met  that  are  20-100  times  drove  physiological  levels 
are  necessary  to  cause  growth  inhibition  in  cultured  cells.  Excessively 
high  concentrations  of  Se-Met  could  produce  a  spectrum  of  nonspecific 
effects  that  may  not  be  related  to  the  anticancer  effect  of  selenium.  To 
obviate  this  problem,  a  stable  monomethylated  selenium  metabolite, 
MSA  (CH3Se02H),  was  developed  specifically  for  in  vitro  studies  (7). 
We  found  that  premalignant  human  breast  cell  lines  were  sensitive  to 
growth  inhibition  and  cell  cycle  block  by  MSA  at  a  concentration  as  low 
as  2.5  /xm  (8).  In  addition,  Jiang  et  al.  (9)  recently  reported  that  MSA 
induced  apoptosis  in  DU- 145  human  prostate  cancer  cells  at  a  concen¬ 
tration  of  5  p.M.  Sinha  et  al.  (10)  also  showed  that  with  mouse  mammary 
tumor  cells,  a  10-min  exposure  to  5  /xm  MSA  was  sufficient  to  cause  a 
change  in  the  expression  of  a  handful  of  genes  as  detected  by  the  Atlas 
mouse  cDNA  expression  array.  The  concentration  of  selenium  used  in  die 
above  studies  is  within  the  physiological  range  of  selenium  in  the  circu¬ 
lation.  As  expected,  MSA  also  has  excellent  anticancer  activity  in  vivo 
(7).  We  are,  therefore,  confident  that  the  information  obtained  with  MSA 
from  cell  culture  studies  would  be  relevant  to  the  action  of  selenium. 

In  this  study,  we  first  examined  the  dose-dependent  effect  of  MSA 
on  the  growth  of  the  PC-3  human  prostate  cancer  cell  line.  We  then 
showed  that  growth  inhibition  by  MSA  was  likely  attributable  to  a 
combined  effect  on  cell  cycle  block  and  apoptosis.  Next  we  used  the 
oligonucleotide  array  technology  to  gain  further  insight  into  the  gene 
expression  changes  that  might  play  a  role  in  the  regulation  of  these 
cellular  events.  Many  potential  selenium-responsive  genes  were  iden¬ 
tified  by  this  method.  These  genes  fell  into  12  clusters  of  distinct 
kinetics  pattern  of  modulation  by  MSA.  The  early  response  genes 
were  grouped  on  the  basis  of  their  known  functions  in  cell  growth  and 
tumorigenesis.  From  the  array  data,  we  were  able  to  develop  an 
integrated  scheme  of  signaling  pathways  that  might  explain  the  action 
of  selenium  in  blocking  cell  cycle  progression. 

MATERIALS  AND  METHODS 

Selenium  Reagents  and  Cell  Line,  MSA  was  synthesized  as  described 
previously  (7),  Se-Met  was  purchased  from  Sigma  (St.  Louis,  MO).  The  PC-3 
human  prostate  cancer  cells  were  obtained  from  American  Type  Culture 
Collection  (Manassas,  VA).  The  cells  were  cultured  in  RPMI  1640  supple- 
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mented  with  10%  feta!  bovine  serum,  100  units/ml  penicillin,  100  fig/ml 
streptomycin,  and  2  mM  glutamine,  and  maintained  in  an  atmosphere  of  5% 
C02  in  a  37°C  humidified  incubator. 

MTT  Cell  Proliferation  Assay.  The  assay,  which  is  based  on  the  conver¬ 
sion  of  the  yellow  tetrazolium  salt  MTT  to  purple  formazan  crystals  by 
metabolieally  active  cells  (1 1),  provides  a  quantitative  determination  of  viable 
cells.  Cells  were  seeded  in  24-well  plates  at  a  density  designed  to  reach 
70-80%  confluency  at  the  time  of  assay.  At  48  h  after  seeding,  cells  were 
treated  with  various  concentrations  of  Se-Met  or  MSA  in  triplicate.  After  24, 
48,  or  72  h  of  treatment,  200  fXi  of  MTT  was  added  to  each  well  of  cells,  and 
the  plate  was  incubated  for  4  h  at  37°C.  The  MTT  crystals  from  both  attached 
and  floating  cells  were  solubilized  in  isopropanol  and  subjected  to  centrifuga¬ 
tion  to  pellet  the  cellular  debris.  Speetrophotometric  absorbance  of  each 
sample  was  measured  at  570  nm  using  a  Spectra  Microplate  Reader  (SLT- 
Labinstruments  Ges.m.b.H.,  Salzburg,  Austria). 

Cell  Cycle  Analysis.  PC-3  cells  were  plated  at  a  density  of  104  cells/cm2  in 
T75  culture  flasks  and  allowed  to  grow  for  48  h  to  reach  70-80%  confluency. 
Synchronization  of  cells  was  achieved  by  starving  in  serum-free  medium  for 
48  h.  Over  85%  of  cells  were  in  G„  phase  at  the  end  of  this  time  period.  On 
returning  to  regular  growth  medium  for  6  h,  cells  were  exposed  to  1 0  ptM  MSA. 
The  procedure  of  serum-starvation  and  refeeding  has  been  described  previ¬ 
ously  by  Sinha  and  Medina  (12)  and  Sinha  el  at.  (13)  to  study  the  effect  of 
selenium  on  cell  cycling.  After  treatment  for  24,  32,  or  48  h,  cells  were 
trypsinized,  washed  in  PBS,  and  fixed  overnight  in  70%  ethanol  at  4°C.  The 
ethanol  solution  was  subsequently  removed  after  centrifugation,  and  cells  were 
resuspended  in  a  buffer  containing  10  niM  Tris  (pH  7.5),  125  mM  sucrose,  2.5 
mM  MgCl2, 0.1 85%  NP40, 0.02  mg/ml  RNase  A,  0.05%  sodium  citrate,  and  25 
jug/ml  PI.  After  incubation  on  ice  for  1  h,  cells  were  analyzed  for  DNA  content 
using  a  FACScan  cytometer  (Becton  Dickinson). 

BrdUrd  Labeling  Assay.  PC-3  cells  were  plated  at  a  density  of  104 
cells'cm2  in  T75  culture  flasks  and  synchronized  as  described  above.  On 
returning  to  regular  growth  medium  for  6  h,  cells  were  exposed  to  10  /x.M  MSA 
for  24  or  for  48  h.  During  the  last  30  min  of  MSA  treatment,  cells  were  labeled 
with  10  jU.M  BrdUrfl  (10  ju.1  of  1  mM  BrdUrd  was  added  to  each  ml  of  culture 
media).  BrdUrd-labeled  cells  were  trypsinized,  fixed,  treated  with  DNase  I, 
and  stained  with  FITC-conjugated  anti-BrdUrd  antibody  using  the  BrdUrd 
Flow  Kit  from  BD  Pharmigen  (San  Diego,  CA).  Stained  cells  were  then 
quantified  by  flow  cytometry,  and  the  data  were  analyzed  with  the  WinList 
software  (Variety  Software  House,  Topsham,  ME). 

Quantitation  of  Apoptosis  by  Flow  Cytometry.  PC-3  cells  were  plated  at 
a  density  of  104  cells/cm2  in  T175  culture  flasks.  At  48  h  after  seeding,  cells  were 
exposed  to  either  5  or  1 0  /um  MSA  for  48  or  72  h.  Adherent  cells  harvested  by  mild 
ttypsinization  were  pooled  together  with  detached  cells.  Cells  were  stained  with 
biotin-conjugated  Annexin  V,  FITC-conjugated  streptavidin,  and  PI  using  the 
Annexin  V-Biotin  Apoptosis  Detection  kit  (Oncogene  Research  Products,  Boston, 
MA)  as  per  the  manufacturer’s  protocol.  Apoptotic  cells  were  subsequently 
counted  by  flow  cytometry,  and  the  data  were  analyzed  with  the  WinList  software 
(Variety  Software  House,  Topsham,  ME). 

Oligonucleotide  Array  Analysis.  PC-3  cells  were  plated  at  a  density  of 
104  cells/em2  in  15-cm  culture  dishes.  Synchronization  was  achieved  as 
described  above.  After  exposure  to  10  p-M  MSA  for  12,  24,  36,  or  48  h,  total 
RNA  and  protein  were  isolated  using  TRIzol  (Life  Technologies,  Inc.).  The 
experiment  was  repeated,  and  the  total  RNA  collected  from  the  replicate  was 
pooled  and  submitted  to  microarray  analysis  using  the  U95A  chip  from 
Affymetrix  (Santa  Clara,  CA).  Biotinylated  eRNA  probe  generation,  as  well  as 
array  hybridization,  washing,  and  staining,  was  carried  out  according  to  the 
standard  Affymetrix  GeneChip  protocol.  Fluorescence  intensity  for  each  chip 
was  captured  with  a  Hewlett-Packard  laser  confocal  scanner.  Absolute  analysis 
of  each  chip  and  comparative  analysis  of  MSA-treated  samples  with  the 
untreated  control  samples  were  performed  by  using  the  Affymetrix  Microarray 
Suite  software.  The  mean  hybridization  signal  for  each  sample  was  set  as  1000 
arbitrary  units  to  normalize  the  signal  values  of  all  of  the  genes  on  the  chip 
(global  normalization)  between  different  samples.  A  treatment/control  signal 
ratio  of  &2  or  S0.5  was  chosen  as  the  criterion  for  induction  or  repression, 
respectively.  These  threshold  values  are  commonly  used  in  the  literature  for 
microarray  expression  analysis  (14-16).  GENECluster  program  (Massachu¬ 
setts  Institute  of  Technology,  Boston,  MA)  and  Affymetrix  Data  Mining  Tool 
were  used  for  clustering  analysis. 


Western  Blot  Analysis.  Western  blot  analysis  was  performed  as  described 
previously  (17)  using  the  TRIzol  isolated  protein.  Briefly,  —50  /xg  of  protein 
was  resolved  over  10-15%  SDS-PAGE  and  transferred  to  polyvinylidine 
difluoride  membrane.  The  blot  was  blocked  in  blocking  buffer  [5%  nonfat  dry 
milk,  10  mM  Tris  (pH  7.5),  10  mM  NaCI,  and  0.1%  Tween  20]  overnight  at 
4°C,  incubated  with  the  primary  antibody  at  37°C  for  1  h,  followed  by 
incubation  with  an  antimouse  antirabbit,  or  antisheep  horseradish  peroxidase- 
conjugated  secondary  antibody  (Bio-Rad,  Hercules,  CA)  at  37°C  for  30  min. 
Individual  protein  bands  were  visualized  by  an  enhanced  chemiluminescence 
kit  obtained  from  Amersham  Pharmacia  Biotech  (Piscataway,  NJ),  Immuno- 
reactive  bands  were  quantitated  by  volume  densitometry  using  the  ImageQuant 
software  (Molecular  Dynamics,  Sunnyvale,  CA),  and  normalized  to  actin.  The 
following  monoclonal  antibodies  were  used  in  this  study  (source):  anti-actin 
(Sigma,  St.  Louis,  MO);  anti-DHFR,  CDK1,  and  CDK2  (BD  Transduction 
Laboratory,  San  Jose,  CA);  anti-PCNA  (Santa  Cruz  Biotechnology,  Santa 
Cruz,  CA);  and  anti-cyclin  A,  cyclin  E2,  CDK4,  p21WAF  1  (NeoMarkers, 
Fremont,  CA).  Polyclonal  antibodies  to  CHK2  and  GADDI 53  were  obtained 
from  Calbiochem  (La  Jolla,  CA)  and  Santa  Cruz  Biotechnology  (Santa  Cruz, 
CA),  respectively. 

Statistical  Analysis.  The  Student’s  two-tailed  t  test  was  used  to  determine 
significant  differences  between  treatment  and  control  values,  and  P  <  0.05  was 
considered  statistically  significant. 

RESULTS 

Sensitivity  of  Human  Prostate  Cancer  Cells  to  MSA.  The  inhib¬ 
itory  effects  of  MSA  and  Se-Met  on  the  accumulation  of  PC-3  cells 
were  assessed  by  the  MTT  assay.  As  shown  in  Table  1,  MSA  was  able 
to  significantly  suppress  the  growth  of  PC-3  cells  in  a  time-  and 
dose-dependent  manner.  At  72  h  of  treatment,  5  fiM  MSA  reduced  cell 
number  by  ~25%.  Increasing  the  concentration  of  MSA  to  10  fiu 
resulted  in  a  more  pronounced  effect,  leading  to  a  greater  magnitude 
of  growth  inhibition  in  a  shorter  period  of  exposure.  In  contrast,  a 
concentration  of  200  or  400  jxm  Se-Met  was  required  to  produce 
significant  decreases  in  cell  number  at  72  h  or  48  h,  respectively.  It  is 
thus  evident  that  MSA  is  much  more  potent  than  Se-Met  in  inhibiting 
growth  of  these  prostate  cells. 

Cell  Cycle  Block  by  MSA.  To  determine  whether  the  decrease  in 
cell  number  accumulation  by  MSA  was  related  to  cell  cycle  arrest,  we 
proceeded  to  assess  the  evidence  of  cell  cycle  perturbation  by  flow 
cytometry  of  ethanol-permeabilized  cells  stained  with  PI.  Synchro¬ 
nized  PC-3  cells  were  treated  with  1 0  /j.m  MSA  for  24,  32,  or  48  h. 
MSA  did  not  cause  any  significant  change  in  cell  cycle  distribution 
(Fig.  1).  However,  flow  cytometry  of  BrdUrd-labeled  cells  showed 
that  MSA  treatment  resulted  in  a  drastic  decrease  in  the  number  of 
cells  synthesizing  DNA  (Fig.  2).  Therefore,  the  data  suggest  that  MSA 
probably  blocked  cell  cycle  progression  at  multiple  stages.  It  should 
be  noted  that  flow  cytometry  of  Pi-stained  cells  would  not  be  able  to 
detect  a  change  in  the  proportion  of  cells  in  different  phases  of  the  cell 


Tabic  1  Effect  of  MSA  or  Se-Met  on  the  accumulation  of  PC-3  cells  at  three 
treatment  durations 


%  of  untreated  control" 

Treatment 

24  h 

48  h 

72  h 

MSA  (/am) 

1 

97.0  ±  3.5 

98.9  +  5.0 

99,0  ±  4.9 

2.5 

107.2  +  5.0 

106.4  ±  3.8 

101,8  +  6.2 

5 

104.2  +  5.0 

103.7  +  5.3 

74.2  +  6.5* 

10 

101.6  +  6.3 

46.3  ±  4.4® 

38.4  +  1.7* 

Se-Met  (/am) 

25 

101.3  +  4.0 

H2.7  ±  1.7 

102.4  +  3,5 

50 

101.5  +  3.3 

110,2  +  2,5 

94.6  +  4.5 

too 

100.1  +  2.1 

110,7  +  1,3 

87.8  +  5.6* 

200 

94.1  +  1.2 

97.2  ±  6.4 

69.3  +  6.4* 

400 

92.8  +  5,5 

75.3  +  4.4* 

50.1  +  3,6* 

a  Results  are  expressed  as  mean  ±  SE  (n  —  4  independent  experiments). 
b  Significantly  different  compared  to  the  corresponding  control  value  (P  <  0.05). 
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Fig.  1.  Cell  cycle  distribution  in  PC-3  cells  treated  with  MSA.  Results  are  expressed 
as  means  ±  SE  (n  =  3).  *,  statistically  significant  ( P  <  0.05)  versus  untreated  control. 


24  hr  48  hr 

Fig.  2,  BrdUrd  ( BrdU )  labeling  of  PC-3  cells  treated  with  MSA.  Results  arc  expressed 
as  means  ±  SE  (n  =  3).  *,  statistically  significant  (P  <  0.05)  versus  untreated  control. 


cycle  ( i.e the  percentage  of  cells  in  G0-G, ,  S  phase,  or  G2-M)  if  there 
is  a  persistent  slowdown  in  the  transition  in  all  phases  of  the  cell  cycle. 

Induction  of  Apoptosis  by  MSA.  In  an  attempt  to  determine 
whether  MSA  might  also  induce  cell  death,  we  incubated  exponen¬ 
tially  growing  PC-3  cells  for  48  or  72  h  in  the  presence  of  5  or  10  /am 
MSA,  and  quantified  the  extent  of  apoptosis  by  flow  cytometric 
analysis  of  cells  labeled  with  annexin  V  and  PI.  Phosphatidylserine 
extemalization  is  a  characteristic  of  cells  undergoing  apoptosis.  An¬ 
nexin  V  has  a  strong  affinity  for  phosphatidylserine.  Staining  cells 
simultaneously  with  annexin  V  and  PI  allows  the  resolution  of  intact 
cells  (double  negative),  early  apoptotic  cells  (annexin  V-positive  and 
Pi-negative),  and  late  apoptotic  or  necrotic  cells  (double  positive), 
which  can  be  located  in  the  lower  left,  lower  right,  and  upper  right 
quadrants  of  the  eytograms  of  Fig.  3 A,  respectively.  Because  only 
cells  that  are  annexin  V-positive  and  Pi-negative  are  truly  represent¬ 
ative  of  apoptotic  cells,  the  percentage  of  this  cell  population  was 
quantitated  from  four  individual  experiments  and  shown  in  a  bar 
graph  form  in  Fig.  35.  MSA  caused  an  induction  of  apoptosis  at  the 
48-h  time  point,  and  the  effect  was  maintained  with  longer  exposure 
to  MSA  for  72  h.  The  increase  of  apoptosis,  which  followed  the 
occurrence  of  growth  arrest,  appeared  to  maximize  with  5  /xm  MSA, 
and  no  further  enhancement  was  detected  with  10  p.M  MSA. 

Profiling  of  MSA-responsive  Genes  by  Oligonucleotide  Array 
Analysis.  We  used  the  Affymetrix  human  genome  U95A  Chip  to 
profile  the  changes  in  gene  expression  and  to  characterize  selenium- 
responsive  targets  that  might  lead  to  cell  growth  inhibition  by  MSA. 
This  GeneChip  contains  probes  to  12,000  known  genes.  As  it  would 
be  more  informative  to  do  the  profiling  at  a  series  of  time  points  than 


to  do  replicate  analysis  at  a  single  time  point,  we  decided  to  commit 
our  available  resources  to  the  former  design.  For  each  time  point  at  1 2, 
24,  36,  and  48  h  post-MSA,  three  separate  preparations  of  RNA 
samples  were  pooled  and  submitted  to  array  hybridization. 

Pairwise  comparative  analysis  between  MSA-treated  samples  and  the 
corresponding  untreated  control  samples  at  each  time  point  was  per¬ 
formed  by  using  the  Affymetrix  Microarray  Suite  software.  This  software 
determines  whether  a  given  gene  is  differently  expressed  based  on  a 
decision  matrix  including  the  net  change  in  intensity  values,  fold  of 
change,  and  other  parameters.  A  no-change  decision  call  was  assigned  a 
value  of  “1.”  Genes  with  expression  changes  of  >2  or  £0.5  were 
considered  as  MSA-responsive  genes.  The  2-fold  difference  limit  was 
chosen  based  on  our  previous  experience  with  microarray  data  analysis 
and  was  also  in  general  agreement  with  other  reported  array  experiments. 
Table  2  shows  the  number  of  genes  induced  or  repressed  by  MSA  at  each 
time  point.  There  were  significantly  fewer  MSA-modulated  genes  at  the 
48-h  time  point  than  at  the  other  three  early  time  points.  This  could  be 
because  growth  inhibition  by  MSA  has  reached  ~50%  at  48  h  (Table  1) 
and  because  the  underpinning  molecular  changes  have  already  peaked 
and  receded  by  this  time. 

To  study  in  detail  the  kinetics  of  expression  changes  in  response  to 
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Fig.  3.  Quantitation  of  apoptotic  cells  by  flow  cytometric  analysis  of  MSA-treated 
PC-3  cells  labeled  with  annexin  V  and  PI.  A,  eytograms  from  flow  cytometric  analysis. 
Intact  cells,  early  apoptotic  cells,  and  late  apoptotic  and  necrotic  cells  are  located  in  the 
lower  left,  lower  right ,  and  upper  right  quadrants  of  the  eytograms,  respectively.  B, 
percentages  of  early  apoptotic  cells.  Data  are  presented  as  means  ±  SE  (n  =  4).  *, 
statistically  significant  (P  <  0.05)  versus  untreated  control. 


Table  2  Number  of  genes  modulated  by  MSA 


Time  point  (h) 

Genes 

12 

24 

36 

48 

Induced  by  MSA 

502 

926 

255 

133 

Repressed  by  MSA 

364 

496 

588 

136 
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Fig,  4.  Average  gene  expression  profiles  for  each  SOM  cluster.  The  fold  of  change  for  each  gene  in  a  cluster  was  averaged  and  plotted  against  duration  of  treatment.  Data  are 
presented  as  means  ±  SE.  Numbers  in  parentheses,  the  number  of  genes  included  in  each  cluster.  Twelve  clusters  were  derived  from  the  clustering  analysts,  demonstrating  different 
kinetics  of  modulation  by  MSA. 


MSA  treatment,  the  genes  modulated  by  MSA  at  one  or  more  time  points 
were  subjected  to  clustering  analysis  using  the  SOM  algorithm.  Cluster¬ 
ing  analysis  avails  to  the  grouping  of  genes  according  to  similarities  in 
their  expression  profiles  across  multiple  time  points.  The  SOM  algorithm 
is  ideal  for  pattern  discovery  and  has  been  found  to  be  useful  in  the 
elucidation  of  biological  pathways,  because  genes  that  are  regulated  in  a 
coordinated  fashion  are  often  related.  A  total  of 2647  genes  were  included 
in  this  analysis.  This  represents  2 1  %  of  the  genes  on  the  U95  A  GeneChip. 
These  genes  fell  into  12  clusters  of  distinct  kinetics  pattern  of  modulation 
by  MSA  (Fig.  4),  with  the  early  response  genes  in  clusters  1, 7,  and  8;  the 
intermediate  response  genes  in  clusters  2,  3,  6,  10,  and  1 1;  and  the  late 
response  genes  in  clusters  4, 5, 9,  and  12.  Because  cell  growth  inhibition 
by  10  /am  MSA  occurred  between  24  to  48  h  (Table  1),  we  consider  the 
genes  modulated  by  MSA  at  the  12-h  time  point  as  the  early-response 
genes.  We  also  believe  that  the  change  in  expression  of  these  early- 
response  genes  is  important  in  initiating  the  cascade  of  events  leading  to 
the  action  of  MSA.  We,  therefore,  classify  these  early-response  genes 
according  to  their  known  functions  in  cell  growth/tumorigenesis 
(Table  3). 

There  is  a  prominent  group  of  genes  implicated  in  cell  cycle 
regulation.  Many  negative  cell  cycle  regulators  were  induced  by 
MSA,  including  checkpoint  proteins  RAD9  and  CHK2,  CDK  inhibi¬ 
tors  pl¥HK4d  and  p21WAe> ,  RB-binding  protein  1,  GADD45,  and 
protein  phosphatase  2C,  On  the  other  hand,  numerous  cyclins,  CDKs, 
and  genes  required  for  DNA  replication  or  mitosis  were  repressed  by 
MSA.  Because  these  genes  control  the  transition  of  different  phases  of 
the  cell  cycle,  the  change  in  their  expression  could  mediate  the 
inhibitory  effect  of  MSA  on  cell  cycle  progression.  There  is  a  second 
group  of  genes  involved  in  apoptosis.  Both  toll-like  receptor  2  and 
caspase  9  were  up-regulated  by  MSA.  Toll-like  receptor  2  is  a  death 
receptor  (18),  and  caspase  9  is  an  upstream  activator  of  the  caspase 
cascade  (19).  In  addition,  MSA  down-reguiated  the  expression  of 
survivin,  a  member  of  the  IAP  (inhibitor  of  apoptosis)  family.  The 
reason  that  so  few  apoptosis-regulatory  genes  were  affected  by  MSA 
at  the  1 2-h  time  point  is  probably  because  the  induction  of  apoptosis 


by  MSA  did  not  become  evident  until  after  48  h.  An  analysis  of 
apoptosis  gene  expression  changes  at  the  24-  and  36-h  time  points  will 
be  the  subject  of  a  separate  study. 

In  Table  3  is  a  group  of  genes  encoding  cytoskeleton  components, 
cell  membrane  glycoproteins,  as  well  as  matrix  metalloproteinases. 
The  responses  of  cytoskeleton  genes  to  MSA  were  varied.  However, 
the  up-regulation  of  invasion  suppressors  (e.g.,  cadherins)  and  the 
down-regulation  of  invasion  activators  (e.g.,  integrins,  endonexin, 
hyaluronan  receptors  CD44  and  RHAMM,  and  MMP21/22)  suggest  a 
possible  role  of  MSA  in  inhibiting  tumor  cell  invasion.  The  table  also 
shows  a  group  of  signal  transduction  genes  that  are  responsive  to 
MSA.  In  particular,  there  is  a  cluster  of  small  GTPases  and  their 
associated  factors,  such  as  Ras-Hke  protein  7c/ 0,  GTPase  activating 
factor-2,  RAN-binding  protein  8,  G  protein-coupled  receptor  37, 
RAB31,  RAB28,  RAB  7 -like  /,  regulator  of  G protein  signaling  10,  Rho 
E,  Rho  2,  and  prenylated  RAB  acceptor  1.  These  genes  belong  to  the 
Ras  and  Rho  family,  the  members  of  which  are  known  to  regulate 
diverse  cellular  functions,  such  as  cell  cycle  progression,  actin  cy¬ 
toskeleton  organization,  malignant  transformation,  and  MAPK  signal¬ 
ing  cascades.  In  addition,  MSA-mediated  up-regulation  of  several 
MAPK  cascade  genes,  including  MEK1,  MEK3b,  MEK5,  and  JNK1, 
may  amplify  its  effect  on  Ras/Rho  signaling.  Another  noteworthy 
observation  is  the  repression  of  a  key  player  of  the  survival  pathway, 
PI3-kinase.  This  suggests  that  selenium  not  only  activates  proapo¬ 
ptosis  signals,  but  may  also  suppress  survival  signals  to  augment  the 
stimulus  to  apoptosis. 

MSA  was  found  to  modulate  a  large  group  of  transcription  factors, 
especially  the  zinc-finger  family  proteins  (ZNFs),  the  mye  proteins 
and  associated  factors,  the  ATF/CREB  proteins  and  their  binding 
proteins,  as  well  as  the  inhibitor  of  DNA  synthesis  (Id)  family  pro¬ 
teins.  Many  of  these  transcription  factors  play  critical  roles  in  the 
regulation  of  cell  cycle  progression,  apoptosis,  and  malignant  trans¬ 
formation.  The  change  in  the  expression  of  these  trans-acting  factors 
could  lead  to  an  altered  transcription  of  a  series  of  other  genes.  To 
wrap  up  the  information  summarized  in  Table  3,  two  growth  factors 
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Tabic  3  Functional  classification  of  MSA  early-response  genes 
A  value  >1  represents  induction;  a  value  <t  represents  repression. 

Gene  Modulation 

Cell  cycle 


RAD9 

3.0 

RB-binding  protein  I 

2.0 

p!9mK4i 

2.0 

GADD45 

2.0 

p21*Af' 

3.0 

PP2C-0 

4.0 

PP2C-a2 

2.0 

CHK2 

3.0 

CDCS 

0.6 

CDC7 -related  kinase 

0.5 

M  phase  phosphoprotein  l 

0.4 

ribonucleotide  reductase  Ml  subunit 

0.4 

STKJ5 

0.5 

Mitosin 

0.5 

thymidylate  synthase 

0.2 

chromatin  assembly  factor-l  pi 50  subunit 

0.5 

CDK1 

0.5 

MCM7 

0.3 

replication  factor  C 

0.4 

replication  protein  A  32-kDa  subunit 

0.5 

dihydrofolate  reductase  (DHFR) 

0.2 

PCTAIRE-1 

0.3 

DNA  pol-e  subunit  B 

0.3 

DNA  pol-a 

0.4 

protein  phosphatase  2A  regulatory  subunit  a 

0.2 

Cyclin  A 

0.3 

Cyclin  E2 

0,3 

MCM6 

0.5 

Ki67a 

0.2 

DNA  primase 

0.4 

Kinesin-like  1 

0.2 

MCM3 

0.4 

Lamin  B2 

0.4 

CDK4 

0.5 

MCM5 

0.2 

PCNA 

0.1 

CDK2 

0.5 

Apoptosis 

toll-like  receptor  2 

2.6 

Caspase  9 

2.0 

Survivin 

0.3 

Angiogenesis 

VEGF-C 

0.4 

Protein  synthesis 

elF-4y 

0.5 

Growth  factor 

bFGF 

0.3 

Wnt?a 

0.3 

Tumor  suppressor/Growth  inhibitor 

BRCA2 

3.2 

DLC-J 

3.1 

TGF- (3 

2.9 

TGF-  J3  type  III  receptor 

3.6 

BMP-4 

4.8 

PTEN 

2.0 

Transcription  factor 

BACH 1 

2.2 

Hox5,4 

5.3 

ZNF345 

2,8 

ZNF217 

3.6 

ZNF165 

3.8 

ZNF267 

2,9 

ZNF75 

2.3 

Ring  ZNF 

2.5 

ZNF274 

3.2 

ZNF278 

0.5 

ZNF  X-linked 

0.3 

Kruppel-like  ZNF 

3.5 

CBP/p3 00-interacting  transcativation 

2.5 

CREBl 

3.0 

ATF5 

0,4 

E74-like  factor  3 

2.4 

SRCl 

2.6 

idl 

0.5 

Id3 

0.2 

Forkhead  box  MIB 

0.4 

N-myc  interactor 

0.5 

Myc-associated  ZNF 

0.2 

b-myb 

0.1 

Jun-B 

0.1 

Table  3  Continued 


Gene 

Modulation 

HNF-3a 

0.5 

TAFinO-a 

0.5 

NFkB  pSO 

0.2 

BAF57 

0.5 

HFH-1IA 

0.4 

GADDI53 

8,0 

DNA  repair 

hPMS2 

6.0 

ERCC1 

3.3 

Signal  transduction 

MEK1 

2.3 

MEK3b 

2.1 

MEKS 

6,0 

JNK1 

3.0 

Ras-like  protein  TclO 

3.8 

GTPase  activating  factor-2 

4.4 

RAN-binding  protein  8 

3.1 

G  protein-coupled  receptor  37 

2.3 

RAB31 

2.5 

RAB2S 

2.1 

RAB7-like  1 

2.4 

Regulator  of  G  protein  signaling  10 

2.0 

Rho  E 

3.2 

Rho  2 

0,5 

prenylated  RAB  acceptor  1 

0.4 

c AMP-regulated  guanine  nucleotide  exchange 

2.2 

factor  II 

CDP-DAG  synthase 

2.3 

phospholipase  C  (34 

2.1 

receptor  of  retinoic  acid 

5.0 

dual  specificity  phosphatase  MKP-5 

6.4 

calmodulin-dependent  protein  kinase  IV 

2.0 

PKC-a 

3.4 

PPI 

2.3 

PP5 

0.2 

PI  3-kinase 

0.4 

Cytoskelcton 

Adducin-y 

2.4 

Calponin 

5.3 

Actin-binding  protein  57 

11.9 

Cytokeratin  20 

4.1 

Tau 

0.2 

Filamin 

0.2 

Tubulin  a  1  isoform  44 

0.4 

Non-muscle  a-actinin 

0.4 

Adhesion/Invasion 

Cadherin-1 5 

6.5 

Cadherin-7 

2.5 

integrin  (3 -5 

0.2 

integrin  (3-4 

0.3 

(3  3-endonexin 

0.5 

CD44 

0.5 

RHAMM 

0.2 

MMP2I/22 

0.1 

(bFGF  and  Writ 7a),  an  angiogenesis  molecule  ( VEGF-Cj ,  and  one 
translation-initiation  factor  gene  ( eIF-4y )  that  is  amplified  in  cancer 
cells  (20),  were  down-regulated  by  MSA.  In  contrast,  three  tumor 
suppressor  genes  ( BRCA2 ,  DLC-I,  and  PTEN),  three  TGF-/3  family 
members  or  receptor  {TGF- j3,  TGF-/3  type  III  receptor,  and  BMP-4), 
and  two  DNA  repair  genes  ( hPMS2  and  ERCC1)  were  up-regulated 
by  MSA.  The  regulation  of  these  genes  may  represent  additional 
mechanisms  by  which  MSA  exerts  its  anticancer  effect. 

Confirmation  of  Array  Data.  We  used  Western  blot  analysis  to 
confirm  the  changes  in  expression  of  a  subset  of  10  cell  cycle  genes; 
CHK2,p21WAF',  GADDI 53,  cyclin  A,  DHFR,  CDK1,  CDK2,  CDK4, 
PCNA,  and  cyclin  E2.  As  shown  in  Fig.  5  and  Table  4,  the  Western 
expression  changes  of  7  genes  {CHK2,p2lWAFt ,  GADDI 53,  cyclin  A, 
DHFR,  CDK1,  and  CDK2)  correlated  well  with  the  array  data.  This 
represents  an  agreement  rate  of  70%.  The  lack  of  complete  concord¬ 
ance  could  be  attributable  to  either  false  positive  signals  of  the  array 
data  or  the  discrepancy  between  transcript  and  protein  expression. 
One  noteworthy  finding  is  that  although  both  array  and  Western 
analyses  showed  a  down-regulation  of  CDK2  by  MSA,  the  Western 
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Fig.  5.  Confirmation  of  array  data  by  Western  blot  analysis.  TRlzol-isolatcd  proteins 
from  PC-3  cells  were  subjected  to  immunoblotting  using  an  antibody  specific  for  CHK2, 
cyclin  A,  p21WAF  l,  GADDI 53,  DHFR,  CDK2,  or  CDK.1.  Signals  were  normalized  to  the 
ones  for  actin  to  control  for  loading  variation.  The  results  shown  here  are  representative 
of  that  from  three  similar  independent  experiments. 


data  additionally  revealed  a  reduction  of  phosphorylated  CDK2  (Fig. 
5;  Table  4).  This  dichotomy  clearly  reinforces  a  fundamental  limita¬ 
tion  of  the  array  technology  in  that  changes  beyond  the  step  of  gene 
transcription  are  not  detectable  by  this  method. 

DISCUSSION 

Prostate  cancer  is  the  most  common  cancer  diagnosed  and  the 
second  leading  cause  of  cancer-related  deaths  in  men  in  the  United 
States.  However,  little  is  known  about  the  etiological  factors  for  this 
disease.  Consequently,  it  is  not  possible  to  institute  primary  interven¬ 
tion  strategies  to  remove  the  causative  agents  from  the  environment. 
Secondary  intervention  strategies  are,  therefore,  necessary  to  reduce 
the  morbidity  and  mortality  of  prostate  cancer,  and  selenium  inter¬ 
vention  has  been  championed  as  a  viable  option.  The  present  study 
was  undertaken  to  investigate  the  molecular  targets  and  the  signaling 
pathways  underlying  the  anticancer  activity  of  selenium  in  prostate. 

Previous  studies  by  Sinha  and  Medina  and  by  Sinha  et  al.  (12,  13) 
showed  that  selenium  is  able  to  block  cell  cycle  progression  at  specific 
checkpoints,  which  might  be  explained  by  a  decrease  in  CDK2  kinase 


activity.  These  experiments  were  done  with  mouse  mammary  tumor 
cells  treated  with  50  /hm  of  methylselenocysteine;  this  concentration 
of  selenium  is  at  least  10  times  higher  than  that  found  in  the  circula¬ 
tion  under  noimal  physiological  condition.  For  this  reason,  MSA  is  a 
more  appropriate  agent  for  in  vitro  studies.  Our  results  demonstrate 
that  MSA  inhibits  the  growth  of  prostate  cancer  cells  by  cell  cycle 
blockade  and  apoptosis.  We  used  the  GeneChip  technology  to  profile 
selenium-mediated  gene  expression  changes  in  a  time  course  experi¬ 
ment.  Of  a  total  of  12,000  genes  screened,  over  2,500  were  identified 
to  be  responsive  to  selenium  treatment.  The  shear  magnitude  of  this 
number  is  somewhat  unexpected.  These  genes  were  grouped  into 
early-,  intermediate-,  and  late-response  clusters.  Because  the  early- 
response  genes  are  likely  to  be  more  important  in  initiating  the  effects 
of  selenium,  we  focused  our  attention  on  them  in  our  follow-up 
analysis. 

Certain  key  cell  cycle  regulators  are  among  the  early-response 
genes.  On  the  basis  of  their  altered  expression,  we  propose  a  number 
of  tentative  signaling  pathways  (in  a  cartoon  format)  that  might 
mediate  the  outcome  of  cell  cycle  blockade  by  selenium.  As  shown  in 
Fig.  6,  selenium  treatment  increases  the  expression  of  p21WAFI,  which 
has  dual  functions  in  regulating  the  activity  of  CDK/cyclin  complexes. 
Although  p21WAFI  is  a  potent  inhibitor  of  cyclin  E/A-dependent 
CDK1/2,  it  promotes  the  assembly  and  the  nuclear  translocation  of 
cyclin  D-CDK4/6  complexes,  leading  to  an  increase  in  cyclin  D- 
associated  kinase  activity  (21).  However,  the  induction  ofp!9INIC4d  by 
selenium  counteracts  the  latter  effect.  The  pl91NK4d  protein  binds  to 
and  inhibits  the  cyclin  D-CDK4/6  complexes,  thus  releasing  p21w  AF1 
from  CDK4  and  CDK6.  The  cooperative  action  of  p  1 9,N  K4  d  and 
p21WAFI  leads  ultimately  to  an  inhibition  of  both  cyclin  D-  and  cyclin 
A/E-dependent  kinases.  The  down-regulation  of  CDK1,  CDK2  and 
cyclin  A  by  selenium  provides  an  amplified  effect  on  this  cascade  of 
events.  Complete  phosphorylation/inactivation  of  pRB  requires  the 
sequential  actions  of  cyclin  D-CDK4/6  and  cyclin  E-CDK2  (22). 
Thus,  p 1 9lN k4 a _  and  p2lWAFI -mediated  inhibition  of  CDK2,  CDK4, 
and  CDK6  could  result  in  decreased  phosphorylation  of  pRB,  Hy- 
pophosphorylated  pRB  interacts  with,  and  negatively  regulates,  the 
activity  of  E2F  transcription  factors.  Loss  of  E2F  activity  prevents  the 
transcription  of  genes,  e.g.,  DHFR  and  cyclin  A,  required  for  progres¬ 
sion  into  S  phase.  In  addition,  although  not  depicted  in  Fig.  6, 
P21WAFI  is  able  to  bind  to  PCNA  and  directly  inhibit  its  activity  (23), 
and  interact  with  E2F  subunits  and  disrupt  E2F-CDK-pl07  DNA 
binding  complex  (24, 25).  These  changes,  collectively,  are  expected  to 
result  in  a  blockade  of  DNA  replication. 

As  shown  in  Fig.  7,  the  elevated  expression  of  CHK2  by  selenium 
treatment  leads  to  increased  phosphorylation  of  CDC25  proteins, 
which  subsequently  bind  to  14-3-3  proteins  and  are  exported  to  the 
cytoplasm.  CDC25  proteins  are  responsible  for  removing  the  inhibi¬ 
tory  phosphates  from  CDK1  and  CDK2,  allowing  them  to  be  activated 


Tabic  4  Comparison  of  expression  changes  detected  by  array  and  Western  analyses 


The  values  represent  treatment:eontrol  ratio. 


Gene 

Array  analysis 

Western  analysis" 

12  h 

24  h 

36  h 

48  h 

DifT  Call* 

12  h 

24  h 

36  h 

48  h 

DifT  Call 

CHK2 

3 

1 

1 

1 

t 

3 

4 

1 

10 

f 

pllWAF‘ 

3 

1 

1 

1 

T 

3 

5 

4 

3 

t 

GADDIS} 

8 

12 

6 

2 

T 

14 

6 

S 

1 

t 

cyclin  A 

0.3 

1 

0.3 

0.4 

1 

0.3 

0.3 

0.4 

1 

1 

DHFR 

0.25 

1 

0.5 

1 

1 

0.3 

0.3 

0.3 

0,2 

1 

CDKI 

0.5 

0.5 

1 

1 

l 

0.3 

0.7 

1 

1 

l 

CDK1 

Phosphorylated 

t 

0.4 

0.06 

0.6 

0.09 

l 

Unphosphorylated 

0.3 

0.3 

10 

0.8 

1 

a  The  value  represents  the  mean  of  three  experiments.  The  SE  is  about  10%. 
h  DifT  Call,  difference  call;  f  ,  increase;  i  ,  decrease. 
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Fig.  6.  Schematic  representation  of  signaling  pathways  down¬ 
stream  of  pl9iNK4j  and  p21WAF1  mediating  MSA-induced  cell 
cycle  blockade  (see  “Discussion").  The  numbers  in  parentheses 
indicate  the  treatmentxontrol-signal  ratio  detected  by  array  anal¬ 
ysis,  Up-regulation  is  denoted  by  a  ratio  of  ^2,  and  down-regu¬ 
lation  is  denoted  by  a  ratio  of  ^0.5.  All  changes,  except  the  one  for 
pl91NK‘,*t>  were  confirmed  by  Western  blot  analysis. 


stimulatory  modification 


transcriptional  modification 


by  CDK  activating  kinases  (26,  27).  A  deficiency  in  nuclear  CDC25 
proteins  prevents  the  activation  of  CDK1  and  CDK2  to  facilitate  cell 
cycle  progression.  PP2Cs  are  phosphatases  known  to  remove  the 
activation  phosphates  from  CDK1  and  CDK2.  An  increased  expres¬ 
sion  of  PP2Cs  by  selenium  will  essentially  lead  to  less  active  CDK1 
and  CDK2.  CDK2  causes  increased  phosphorylation  of  MCMs  either 
directly  or  indirectly  through  CDC7.  MCMs  are  important  for  transi¬ 
tion  to  S  phase.  By  decreasing  the  expression  of  CDK2,  CDC7,  and 
MCMs,  selenium  is  able  to  block  DNA  synthesis.  Selenium  can  also 
induce  a  G2-M  block  by  increasing  the  expression  of  the  checkpoint 
protein,  RAD9.  Furthermore,  the  expression  of  GADD153  is  known 
to  increase  in  response  to  a  variety  of  growth  arrest  or  DNA  damage 
signals  (28-35).  GADDI 53  plays  an  essential  role  in  cell  cycle 


control  and  apoptosis  (36).  We  recently  found  that  in  premalignant 
human  breast  cells,  selenium  induced  an  8-fold  increase  in  the  ex¬ 
pression  of  the  GADDI 53  gene  (8).  Similarly,  a  6-14-fold  induction 
(Table  4)  was  also  observed  in  PC-3  cells.  In  summary,  Figs.  6  and  7 
show  that  different  pathways  modulated  by  selenium  all  converge  to 
block  cell  cycle  progression. 

It  is  clear  that  selenium  affects  not  just  one  key  target,  but  a 
multitude  of  targets.  In  doing  so,  the  impact  of  selenium  is  amplified. 
The  diversity  of  the  molecular  responses  also  makes  it  difficult  for 
transformed  cells  to  escape  the  inhibitory  effect  of  selenium.  A 
reassuring  aspect  of  our  results  is  the  considerable  overlap  of  the 
selenium-modulated  genes  or  signaling  pathways  identified  in  pros¬ 
tate  cells  with  those  previously  identified  in  breast  cells  (8).  The 


Fig.  7.  Schematic  representation  of  signaling 
pathways  downstream  of  CHK2,  GADDI 53,  and 
RAD9  mediating  MSA-induced  cell  cycle  blockade 
(see  “Discussion").  The  numbers  in  parentheses 
indicate  the  trcatmentxontrol-signal  ratio  detected 
by  array  analysis.  Up-rcgulation  is  denoted  by  a 
ratio  of  S2,  and  down-regulation  is  denoted  by  a 
ratio  of  ^0.5. 
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congruency,  however,  is  highlighted  against  a  backdrop  of  certain 
differences  in  cell  genotype  and  methodological  issues.  The  PC-3 
human  prostate  cancer  cells  are  null  forp53,  whereas  the  premalignant 
MCF1QAT  human  breast  cells  have  a  functional  p53,  suggesting  that 
p53  is  not  required  for  the  action  of  selenium  The  12,000-gene 
oligonucleotide  GeneChip  was  used  in  the  prostate  cell  study,  whereas 
the  200-gene  membrane-based  cDNA  array  was  used  in  the  breast  cell 
study.  Despite  these  differences,  similar  selenium  targets  were  iden¬ 
tified,  including  GADDI 53, cyclin  A,  CDK1,  CDK2,  CDK4,  CDC25, 
E2Fs,  as  well  as  the  MAPK/JNK  and  phosphoinositide-3  kinase 
pathways,  thus  lending  confidence  to  the  array  data.  Previous  studies 
of  gene  expression  changes  in  response  to  selenium  have  generally 
been  limited  to  the  analysis  of  a  few  or  a  subset  of  genes  and, 
therefore,  have  provided  only  a  very  narrow  view  of  the  entire 
landscape.  As  far  as  we  are  aware  of,  this  is  the  most  comprehensive 
gene  expression  profiling  study  on  the  molecular  mechanism  of  sele¬ 
nium  in  chemoprevention. 

Venkateswaran  et  al.  (37)  recently  reported  that  PC-3  cells  are  not 
growth  inhibited  by  Se-Met  unless  they  are  transfected  with  a  func¬ 
tional  androgen  receptor.  Suffice  it  to  point  out  that  the  sensitivity  of 
PC-3  and  another  androgen-independent  prostate  cancer  cell  line, 
DU-145,  to  Se-Met  has  been  documented  by  two  other  groups  of 
investigators  (38,  39).  Furthermore,  Jiang  et  al.  (9)  showed  that 
DU- 145  cells  can  be  induced  to  undergo  apoptosis  by  MSA  at  phys¬ 
iological  concentrations.  Together  with  our  study,  the  weight  of 
evidence  seems  to  favor  the  notion  that  the  responsiveness  of  prostate 
cancer  cells  to  selenium  is  not  dependent  on  the  presence  of  a 
functional  androgen  receptor.  Although  androgen  plays  a  critical  role 
in  prostate  carcinogenesis,  a  significant  proportion  of  prostate  cancers 
eventually  become  androgen-unresponsive  and  refractory  to  hormonal 
therapy.  The  fact  that  androgen-unresponsive  cells  are  sensitive  to 
selenium-induced  growth  inhibition  is  encouraging  because  it  sug¬ 
gests  that  selenium  intervention  may  be  a  viable  strategy  for  prevent¬ 
ing  prostate  cancer  recurrence  after  prostatectomy. 
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